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Abstract 

We present the results of the spectroscopic, photometric and orbital pe¬ 
riod variation analyses of the detached eclipsing binary V482 Per. We derived 
the absolute parameters of the system (Mi = 1.51 Mq, M 2 = 1.29 Mq, Ri 
= 2.39 Rq, R 2 = 1.45 R 0 , Li = 10.15 Lq, L 2 = 3.01 Lq) for the first time 
in literature, based on an analysis of our own photometric and spectroscopic 
observations. We conhrm the nature of the variations observed in the sys¬ 
tem’s orbital period, suggested to be periodic by earlier works. A light time 
effect due to a physically bound, star-sized companion (M 3 = 2.14 Mq) on a 
highly eccentric (e = 0.83) orbit, seems to be the most likely cause. We argue 
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that the companion can not be a single star but another binary instead. We 
calculated the evolutionary states of the system’s components, and we found 
that the primary is slightly evolving after the Main Sequence, while the less 
massive secondary lies well inside it. 
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1. Introduction 


The mass of a star is the single most important parameter which deter¬ 
mines the way it lives and dies and what remains from it after its death. It 
can only be directly and precisely determined with the analysis of the light 
and radial velocity variations occurring due to the orbital motion of the com¬ 
ponents of an eclipsing binary system. Since each component of a binary is 
conhned to the space that we call the Roche lobe, its evolution will be dif¬ 
ferent from that of a single star once it hlls this space. Mass transfer occurs 
between the two after contact has been made and disturbs the evolution of 
the components as single stars. Detached binaries are the systems whose 
components have not hlled their Roche lobes yet and thus evolve as single 
stars. Therefore the components of detached eclipsing binaries with well de¬ 
termined parameters have been regularly used to test the stellar evolution 
models. 

V482 Per (BD-|-47° 961, GSC 3332-314) is a detached eclip sing binary 


which was hrst observed and classified as a short period variable by 


Hoffmeister 


(119661) . iHarvig fc LeisI (119811 ) investigated this star based on archival photo- 
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graphic plates. They classified it as an eclipsing sys tem, fonnd its magnitnde 


and p nblished the first light-cnrve and its elements. 


Agerer &: Lichtenknecker 


fll99ll) . based on their first photoelectric observatio ns, pnb li shed B and V 


light cnrves and the B-V color cnrve of V482 Per. 


Popperl fll996l) fonnd a 


mean spectral type of F2 for V482 Per by nsing eqnivalent widths of the Na- 
D line of both components. In the s ame stndy he pointed to t h e disc repancy 
between the photometric stndy bv lAgerer fc Lichtenkneckeii fjl99ll) . which 
implied a middle-F spectral type for the primary, and the weakness of the 
Na-D lines in his spectra from which he determined an earlier spectral type. 
V482 Per was also fonnd to be a close binary in a triple system from the 
light-time effect o bser ved on its 0-C diagra m. This snbject was discnssed by 


Wolf et al. 


fenndh and 


Qgloza et al 


fj 2012 |) . both of whom gave very similar 


valnes of the basic parameters of the third body (P 3 and e). 

In 2011, we acqnired precise photometric measnrements of the system 
in Gerostathoponlion Observatory of the University of Athens. We also ob¬ 
tained spectroscopic data for the system in order to determine its orbital 
parameters and the ratio of the masses of its components. We p e rformed 


an analysis nsing the Wils on-Devinney ( Wilson &: Devinne 


19901: 


Kallrath et al. 


I 


1971 


Wilson! . 


19981) model of onr light cnrves based on onr spectro¬ 


scopically determined mass ratio of the system. As a resnlt, for the first 
time we derived its absolnte parameters (masses, radii, and Inminosities) 
for both of the spectral t ypes (F2 a nd FS) s nggest ed by earlier works of 


Agerer &: Lichtenknecker! f!l99l! ) and !Popper! fjl9961) . For the best fit we 


needed to inclnde a third light, which wou l d also sup port the triple - system 


hypothesis of earlier works by 


Wolf et al. 


(2004 

) and 

Ogloza et al. 

( 2012 ) 
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We have performed an 0-C analysis to ensure that an expected orbital vari¬ 
ation is actually observed. Consequently, we conhrmed the existence of a 
physically bound, stellar size third body and determined its parameters. 


2. Observations and Data Reduction 


Between October and December 2011 (in 13 nights), we obtained mul¬ 
ticolor photometric data for V482 Per with the 40 cm Cassegrain telescope 
located at the Gerostathopoulion Observatory of the University of Athens. 
The telescope’s focal ratio was converted to f/5.1 from f/8 by using a focal 
reducer. It was equipped with an SBIG ST-10 XME GGD detector and a 
set of Johnson-Gousins BVRI hlters. We present a log of these observations 


in Table-dJ We have corrected all t 


using the G-munipack flHrochl . 


l e GGD ima ges for instrumental effects. 


1998 


Moti 


2004]) software package. Next we 


performed differential apperture photometry with the same software, which 
allowed us to experiment with different aperture sizes until we achieved a 
relatively low scatter. Individual differential magnitudes were derived with 
respect to a chosen comparison star (GSG 3332-2173), the stability of which 
was tested against a check star (GSG 3332-1993). We constructed the light 
curves using these differential magnitudes in each hlter and the light ele¬ 
ments that we calculated from a linear £t to the previously published times 
of minima. 

One of the authors of this study (RHN), carried out spectroscopic ob¬ 
servations at the Dominion Astrophysical Observatory (DAO) in Victoria, 
British Golumbia, Ganada, using the Gassegrain spectrograph attached on 
the 1.85 m Plaskett telescope with a grating (^21181) of 1800 lines/mm. 
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Table 1: Information About The Observations 


Log of The Observations 

Dates of Observation (in 2011) 


17,19,20,31 Oct 


1, 

2, 3, 8, 16, 19, 26 Nov 



1, 4 Dec 

# of Data Points 

1630 (B), 

, 1597 (V), 1614 (R), 1549 (I) 

Mean Errors (c) 

0.0031 (B), C 

).0031 (V), 0.0031 (R), 0.0027 (I) 

Star 

mv (mag) 

B-V (mag) 

Variable (V482 Per) 

10.30 

0.40 

Comparison (GSC 3332-2173) 

9.63 

0.32 

Check (GSC 3332-1993) 

11.1 

0.40 

Light elements employed to phase the observations 

Reference Epoch (in HJD) 

Orbital Period 


2455868.4005 

2*^.44677 


blazed at 5000 A giving a reciprocal linear dispersion of 10 A/mm in the 
first order and covering a wavelength region from 5000 to 5260 A. The re- 
dnctions (cosmic hit removal, median backgronnd htting and snbtraction for 
each wavelength, apertnre snmmation, and continnnm normalization) were 


performed with the ’RaVeRe’ software flNelsonl . 


20101) . The same software 


was nsed for wavelength calibration and linearization nsing the Fe-Ar spectra 
as wavelength standards. The hve spectra we nsed are depicted in Fig. [T] 

In order to obtain radial ve l ocity measnrements we applied the Rn cinski 


broadening fnnctions (IRncinski 


2004 : 


Nelson et ah 


2006 


Nelson! . I 2 OI 0 I) (Fig. 


[2D. We phased these observations with the same light elements that we used 
to calculate orbital phases of our photometric measurements, after having 
converted all the observation times to heliocentric values. We found the 
gamma velocity (V..y) of the system to be -22.72 (±3.05) km/s. The resultant 
radial velocity curves for each of the components and the best hts to them are 
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Figure 1: The spectra used in radial velocity analysis. Four individual spectra are shifted 
arbitrarily for clarity. 
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Figure 2: Broadening functions used to determine radial velocities. Some of the functions 
are shifted arbitrarily for clarity. 

given in Fig. El We determined semi-amplitudes of each of the htted curves 
(Ki and K 2 ) and computed the spectroscopic mass ratio of the system (q = 
m^/mp = K 1 /K 2 = 0.856). 


3. The Light Curve Analysis 


We used t 


Wilsonl. 


19901: 


le Wilson-Dev i nney 


Kallrath et af 


code (WD) flWilson &: Devinnevl. 


1971 


19981) to derive parameters of components of 


V482 Per. Instead of using all the individual observations, approximately 
150 mean points were used to speed up computations. They were calculated 
in such a way that they evenly covered the observed light curve. In order 
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Figure 3: Radial velocity data for V482 Per and the best fits in solid lines. 
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to make sure that the global minimum was found within the set ranges of 
free parameters, we made use of the Monte Carlo search algorithm. We fol¬ 


lowed the procedure outlined in 


Zola et ah! fl2014j) . in which the mass ratio 


is hxed to the spectroscopically determined value. The albedo and gravity 
darkening coefficients were set at their theoretical values for a convective 
envelo pe fO.5 and 0 .32, r espectively) suitable for both F2 and F5 spectral 


types (Ivan Hammel. 


the ta bles by 


1993h. The limb darkening coefficients wer e taken from 


Claret fc Bloemen 


km) 

Claret et ah 

(2012ll aud 

Claret et ah 


(120131) and were built into the code. They were selected according to the 
temperature and the wavelength of the observations. We adjusted the in¬ 
clination, temperature of the secondary, surface potentials, and the lumi¬ 
nosity of the primary component. We mad e check s for a ll three existing 
spectral types of this sys t ems: F2 given by IPopperl (119961) . F5 derived by 


Agerer fc Lichtenknecked (jl99ll) from the B-V color and AO listed in the 


SIMBAD database. Due to a possibility that this system has a compan¬ 
ion, third light was added into the model. We hxed the temperature of the 
primary at the values of 6700 K, 6460 K and 9420 K, cor responding to the 


three spectral types according to the calibration given by 


Harmaned (119881) 


and achieved three solutions. It soon turned out that the model with the 
assumed AO spectral type of the primary resulted in highly overluminous 
values for their masses. Therefore, this solution was discarded. 

The conhguration of the system was found to be detached with both 
components to be well within their Roche lobes. We obtained a good £t 
to the observed light curves with just small discrepances in the phase range 
between 0.30-0.45. We derived a small contribution of the third light, ranging 
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Table 2: Best-fit parameters from Wilson-Devinney light curve modelling. 




F2 


F5 

Stellar Parameters 

Value 

Error 

Value 

Error 

Ti [K] 

6700 

Fixed 

6460 

Fixed 

T2 [K] 

6340 

18 

6130 

17 


5.38 

0.04 

5.39 

0.04 

£^2 

7.42 

0.15 

7.38 

0.15 

q = m, 2 lm\ 

0.856 

Fixed 

0.856 

Fixed 

i [ 0 ] 

83.2 

0.3 

83.2 

0.3 

Luminosities 

Li [B] 

9.68 

0.33 

9.63 

0.34 

Li [V] 

9.07 

0.33 

9.06 

0.34 

Li [R] 

8.92 

0.32 

8.91 

0.33 

Li [I] 

8.79 

0.32 

8.72 

0.34 

L 2 [B] 

2.75 

0.09 

2.80 

0.09 

L 2 [V] 

2.71 

0.08 

2.71 

0.08 

L 2 [R] 

2.74 

0.08 

2.75 

0.08 

L2 [I] 

2.78 

0.09 

2.82 

0.09 

^3 [B] 

0.01 

0.02 

0.01 

0.03 

^3 [V] 

0.05 

0.03 

0.05 

0.03 

h [R] 

0.06 

0.02 

0.06 

0.03 

h [I] 

0.07 

0.02 

0.07 

0.03 


Note: Ti ,2 - temperature of the primary and secondary, Hi ,2 - dimensionless surface potentials of the 
components, q = m^lvrix - the system mass ratio, i - orbit inclination (in degrees), Z/ 1^2 - WD code 
luminosities, Z 3 - third light 
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Figure 4: Comparison between observations of V482 Per (BVRI filters - shown from top 
to bottom - shifted arbitrarily for clarity by the amounts given in the figure) and the best 
fits (solid lines) obtained for the F2 spectral type. 

from 1 to 7%, depending on the passband. 

The results derived from the light curve modeling are presented in Table 
121 while model light curves along with observations are shown in Figs. |1]& 
|5]for both F2 and F5 spectral types, respectively. 


4. Eclipse Timing Analysis 


V482 Per has been known to display orbital peri o d var i ations, which 


were suggested to be of a periodic nature flWolf et al 


20041: 


Qgloza et al. 
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Figure 5: The same as the Fig0]but for the spectral type F5. 


12 



Table 3: Results of the eclipse timing analysis. 


Parameter 

Value 

Error 

P 3 [years] 

16.56 

0.11 

A [days] 

0.018 

0.001 


206 

2 

e 

0.83 

0.02 

a 3 sin{i) [AU] 

4.74 

0.22 

/(ms) [Mq] 

0.39 

0.06 

m 3 [Mq] 

2.14 

0.05 


20121) . We have collected all available eclipse times publis hed in the literature 
( Appendix A ) including these determined from our data ihiakos fc Niarch^ . 
20 111 ) and constructed an 0-C diagram (Fig. [6]). We conhrmed that the vari¬ 
ations are in fact periodic and most likely they can be caused by the light 
time effect due to a third body gravitationally bo und to the system on an 


eccentric orbit. We made use of Irwin’s formalism fllrwinl. 


1952 


19591) during 


our analysis, and obtained the parameters of the third body presented in 
Table-Elas a result. 


5. Results Discussion 

We have analyzed light and radial velocity curves, and the long term or¬ 
bital period variations of the detached system V482 Per. We have derived two 
sets of absolute parameters (Table-Hj) as a result of the hrst ever combined 
analysis of its photometric and spectroscopic obse rvations. One set of solu¬ 


tions was based on the spectral type F5 given by lAgerer fc Lichtenknecker 
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Figure 6: Period variation behavior of V482 Per. The best fit is shown as solid line. 
Residuals from the fit are given in the lower panel. The symbols used have been defined in 
the legends. Times of minima are listed in the [Appendix Altogether with their references. 
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(1991 

) and the other was based on F2 spectral type given by 

PoDoer 

(1996) 


We u sed the web interface to query Geneva Stellar Models’ database flEkstrom et al 


2 OI 2 I) and computed evolutionary tracks for three different masses (1.3, 1.5, 


1.7 Mq) in close proximity to the derived masses of V482 Per’s components 
(Mp = 1.51 Mq, M^ = 1.29). The stellar evolutionary models have been 
taken from a den se grid for low r nass s tars assuming no rotation and for the 
solar metallicity (iMowlavi et al.l. l2012l) . 

We next constructed the Hertzsprung-Russell (HRD), Mass-Luminosity (MLD), 
and Mass-Radius (MRD) Diagrams (Figs. [7] & |8]) based on the parameters 
we obtained from our analysis for V482 Per, and the tracks from the models. 
We have given the positions of both components for each of the two solutions 
in the HRD (Fig. [7]) with hlled (primaries) and unhlled (secondaries) circles 
(F2) and squares (F5). The evolutionary model tracks for the three chosen 
masses are also drawn by solid lines in the HRD. Since masses of the compo¬ 
nents do not depend on the assumed temperature of the primary, and both 
solutions resulted in almost the same values for the radii of components, only 
parameters resulting from the F2 spectral type model are presented in the 
MRD (right panel of Fig |8]). However, due to the different temperatures of 
the components between the two solutions, there is a difference in luminosi¬ 
ties of about 12%, which is shown in the MLD by giving both of the solutions 
(left panel of Fig. [8]). 

The positions of V482 Per’s components in all evolution diagrams in Figs. 

[7] & [8] indicate the primary to be somewhat evolved (just above the TAMS) 
and more luminous than expected for its mass. This statement is true for 
solutions obtained for both spectral types assumed. Based on the position 
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Figure 7: A Hertzsprung-Russell Diagram with the posi tions of selected det ached (pri¬ 
maries in filled, secondaries in unfilled grey circles) from Ixorres et all ( 201(1 '!. Positions 
of the components of the system V482 Per are marked with larger symbols (circles show 
F2-spectral type solution, squares that of F5). The evolutionary tracks for a dense grid 
of mass es, computed by inte rpolation between the existing tracks of the Geneva stellar 


models (Mowlavi et al. 


2012 1 are plotted with solid lines while ZAMS and TAMS are in 


dashed lines. 
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Figure 8: Mass-Luminosity (left panel) and Mass-R adius (right panel) dia grams. ZAMS 


and TAMS computed from the theoretical models by 


Mowlavi et al. 


solid lines. Positions of components in selected detached systems from 
are also shown. Symbols mean the same as in Fig. [T] 


(l2oi^ are plotted by 


Torres et al 


(l2mnh 


of the secondary (see Figs. [7] and [8]), it is evolving on the Main Sequence, 
about half way between ZAMS and TAMS. It seems that a better agreement 
with the evolutionary models was obtained for the hotter, F2 type primary 
model. Within the F5 solution the secondary is signihcantly underlumious 
(see Figs. [7] and [8]). 

Our analysis of the system’s 0-C diagram provides a stro ng evidence 


for th e exis tence of a thi r d bod y, proposed in earlier works of 


Wolf et al. 


(120041) ■ and lOgloza et al.l (120121) . Knowing the absolute parameters of the 
system’s components, and assuming a coplanar orbit with the binary, we 
found the companion to have a very close lower mass limit (2.14 Mq) to their 
estimates based on guessed absolute parameters of the system. The orbital 
period (~16 years), eccentricity (~0.83) and the size of the gravitationally 
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Table 4: Absolute Parameters 


Parameter 

F 2 model 

Value Error 

F5 model 

Value Error 

a[RQ] 

10.76 

0.16 

10.76 

0.16 

Mi[Mo] 

1.51 

0.05 

1.51 

0.05 

M2[Mo] 

1.29 

0.05 

1.29 

0.05 

Ri [Rq] 

2.39 

0.04 

2.40 

0.04 

i?2 [7?©] 

1.45 

0.04 

1.46 

0.04 

Li[Lq] 

10.15 

0.69 

8.85 

0.62 

L2[Le] 

3.01 

0.17 

2.65 

0.15 


Note: a - orbital semi-major axis in units of solar radius, Mi ^2 - masses, Ri ^2 - 
radii, Li ^2 - luminosities of the primary and secondary, respectively. 


bound component’s orbit are also very similar in both of these studies and 
ours. We conhrm their hndings based on the absolute parameters determined 
from combined spectroscopy and photometry analysis. The 2.14 Mq tertiary 
companion would be an Al type star with about 9100 K temperature and a 
radius of 2.02 Rq. Its luminosity would be about 25.5 Lq, twice exceeding 
that of the binary system. However, we found a very small contribution 
reaching about 7% in the I hlter and decreasing towards shorter wavelengths. 
From a crude estimate it follows that the companions should have about 1 Lq 
in the near infrared and smaller in shorter wavelengths. Unfortunately, we 
haven’t seen a trace of a third body in the spectra which might be due to 
our resolution limit. As an explanation, a binary system consising of a dwarf 
and a compact object would agree with our hndings. Further radial velocity 
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observations are needed to verify this hypothesis. 
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Appendix A 


Table A.l: Eclipse Timings for V482 Per 


HJD-2400000 

Type 

Method 

Ref 

HJD-2400000 

Type 

Method 

Ref 

47565.3737 

II 

pe 

1 

53318.9165 

I 

ccd 

18 

47823.5048 

I 

pe 

1 

53656.5685 

I 

ccd 

19 

47840.6360 

I 

pe 

1 

53668.8015 

I 

ccd 

18 

47850.4210 

I 

pe 

1 

53685.9293 

I 

ccd 

18 

47943.4012 

I 

pe 

1 

53739.7567 

I 

ccd 

20 

48222.3268 

I 

pe 

1 

53744.6506 

I 

ccd 

20 

48222.3276 

I 

pe 

2 

53766.6715 

I 

ccd 

20 

48299.4004 

II 

pe 

2 

53793.5857 

I 

ccd 

20 

48606.4669 

I 

pe 

3 

54055.3865 

I 

ccd 

21 

48650.5058 

I 

pe 

3 

54057.8341 

I 

ccd 

20 

49060.3247 

II 

pe 

4 

54073.7380 

II 

ccd 

20 

49372.2657 

I 

pe 

5 

54193.6316 

II 

ccd 

22 

49625.5031 

II 

pe 

6 

54400.3747 

I 

pe 

23 

49761.3046 

I 

pe 

7 

54402.8237 

I 

ccd 

22 

50047.5745 

I 

ccd 

7 

54408.9423 

II 

ccd 

22 

50079.3873 

I 

ccd 

7 

54433.4101 

II 

pe 

23 

50106.2985 

I 

ccd 

8 

54515.3736 

I 

ccd 

24 

50380.3383 

I 

pe 

9 

54764.9392 

I 

ccd 

22 

51185.3277 

I 

pe 

10 

54764.9399 

I 

ccd 

22 

51373.7390 

I 

ccd 

11 

54774.7265 

I 

ccd 

22 

52250.9000 

II 

ccd 

12 

54786.9621 

I 

ccd 

22 

52266.8056 

I 

ccd 

12 

54801.6414 

I 

ccd 

22 

52276.5957 

I 

ccd 

12 

54812.6523 

II 

ccd 

22 

52287.6027 

II 

ccd 

12 

54816.3222 

I 

ccd 

25 

52288.8255 

I 

ccd 

12 

54840.7874 

I 

ccd 

22 

52589.7781 

I 

ccd 

12 

54840.7884 

I 

ccd 

22 

52644.8308 

II 

ccd 

13 

54845.6802 

I 

ccd 

22 

52724.3527 

I 

ccd 

14 

54845.6826 

I 

ccd 

26 

52907.8584 

I 

ccd 

13 

54867.7018 

I 

ccd 

22 

52949.4538 

I 

ccd 

15 

54867.7024 

I 

ccd 

22 

52949.4565 

I 

ccd 

16 

55158.8543 

I 

ccd 

27 

53032.6433 

I 

pe 

17 

55169.8635 

II 

ccd 

27 

53317.6900 

II 

ccd 

18 

55185.7682 

I 

ccd 

27 


24 







Table A.l Continued from the previous page 


55201.6755 

II 

ccd 

27 

55563.7765 

II 

ccd 

27 

55245.7097 

II 

ccd 

27 

55573.5628 

II 

ccd 

27 

55432.8751 

I 

ccd 

27 

55579.6804 

I 

ccd 

29 

55443.8859 

II 

ccd 

27 

55590.6906 

II 

ccd 

27 

55459.7893 

I 

ccd 

27 

55599.2520 

I 

ccd 

30 

55476.9162 

I 

ccd 

27 

55847.6067 

II 

ccd 

31 

55476.9168 

I 

ccd 

28 

55848.8258 

I 

ccd 

32 

55481.8093 

I 

ccd 

27 

55848.8265 

I 

ccd 

32 

55497.7146 

II 

ccd 

27 

55852.4969 

II 

ccd 

31 

55498.9346 

I 

ccd 

27 

55865.9522 

I 

ccd 

32 

55498.9364 

I 

ccd 

27 

55865.9527 

I 

ccd 

32 

55503.8308 

I 

ccd 

27 

55868.3970 

I 

ccd 

31 

55508.7236 

I 

ccd 

27 

55875.7404 

I 

ccd 

32 

55509.9476 

II 

ccd 

27 

55940.5796 

II 

ccd 

32 

55514.8421 

II 

ccd 

27 

55946.6960 

I 

ccd 

32 

55519.7355 

II 

ccd 

27 

56221.9605 

II 

ccd 

33 

55536.8653 

II 

ccd 

27 

56232.9677 

I 

ccd 

34 

55557.6609 

I 

ccd 

27 

56324.7231 

II 

ccd 

33 

55563.7761 

II 

ccd 

27 
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